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We introduce the chemical potential in a system of massless fermions in a bag by impossing 
boundary conditions in the Euclidean time direction. We express the fermionic mean number in 
terms of a functional trace involving the Green's function of the boundary value problem, which we 
, . study analytically. Numerical evaluations are made, and an application to a simple hadron model 

' is discussed. 
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> ■ I. INTRODUCTION 

(N ■ 

o> : 

Functional determinants of elliptic differential operators have wide application in describing one-loop effects in 
OO ' Quantum and Statistical Physics. When a bounded manifold in the Euclidean space-time is considered, boundary 

' conditions have a relevant influence on the behavior of such determinants. 
^ : In previous papers, we studied the relationship between functional determinants of elliptic boundary value problems 
' and the corresponding Green's functions ||l|-^. This allowed for the description of confined fields at finite temperature. 
, In the framework of this functional techniques, we studied, in a simplified scenario in two dimensions, the Gibbs 
' free energy of confined massless fermions Q. The chemical potential fj, was introduced through suitable boundary 
CIh' conditions impossed on the Dirac operator in the Euclidean time direction. 

1^ The present paper is an extension of that work: Here we consider a similar calculation for a four dimensional M.l.T. 

bag model. In this effective model of hadrons, a massless fermionic field is confined into a static sphere of radius R and 
subject to local spatial boundary conditions. This problem has been considered previously j|] through an asymptotic 
expansion of the energy density, for large values of the bag radius. As we will show, the knowledge of the complete 
^ , Green's function for this boundary value problem allows us to describe the mean fermionic number and the Gibbs 
' free energy even at small radius, thus improving previous approaches. 

The behavior of hadronic matter at finite temperature and density has received considerable attention (see also 
and references therein) during the last years. The main motivation behind this effort is an attempt to understand 
properly the generally accepted possibility of a deconfining phase transition from hadronic matter to the quark-gluon 
plasma with applications to relativistic heavy ion collisions and to the early Universe. 

Different techniques have been used in connection with this problem. Among others, the lattice approach the 
use of effective potentials for composite operators chiral perturbation theory |l^, the QCD Sum Rule method 
(see for example p2|), and models like the Skyrmion approach and mixed constructions of hybrid models p"^-p^ can 
be mentioned. 

In section^ we obtain the fermionic mean number in terms of a trace involving the Green's function of the problem. 
This last has a nonregular behavior, which we study analytically in order to extract physical results. 

In section |l| the difference between Gibbs and Helmholtz (u — case) free energy is obtained. Complementing 
this result with the Helmholtz energy previously evaluated in , we discuss its application to a simple quark-gluon 
droplet model P,pT[. 

II. THE GIBBS FREE ENERGY AND MEAN FERMIONIC NUMBER 

In this section we will be interested in the study of the Gibbs free energy of a massless fermionic field confined 
inside a static sphere of radius R and obeying local boundary conditions corresponding to the M. 1. T. bag model. 



X 



1 



This amounts to study the functional determinant of the Dirac operator under such boundary conditions, which is 
related to the Grand canonical partition function according to 

E{T, R, fi) = e-^^(^'«-'^) ~ Det {D{P, R) - i^7"),^ . (1) 

Here 

D{P, R) = ^-f'^dt + ■ V, for < t, \x\ < 1, (2) 

and "6c" means that the operator is defined on a space of functions satisfying the local boundary conditions 

Bi;{t, x) ^{U + i i) ^{t, x)^0, for \x\ = 1, (3) 
at the spatial edge and antiperiodic boundary conditions in the Euclidean time direction, 

=-^(0,f) (4) 

as corresponds to the fermionic character of the field. 

In the above expression, stands for the chemical potential and /3 denotes the inverse of the temperature. In our 
conventions, the Euclidean gamma matrices are 

70 = «> I2 , 7 = ® , 75 = (g) I2 , (5) 

where p*^, cr'^, fc = 1, 2, 3, denote the Pauli matrices. 

Following Q], we can relate G(T, i?, p) to the Green's function of this boundary value problem, 

-/3f (/3, R, /i) = Tr{^ In {D{f3, R) - 

(6) 

= Tr{-ij°Kb,{t,x;t',x')}, 

where Ki,c{t, x; t' , x') satisfies 

{D{I3, R) - ip7°) Kbcit, X- 1', x') = S{x - x')S{t - t') 

BKbait, x; t', x') = 0, for |f | = 1, (7) 
Kbc{l,x-t',x') ^ -Kbc{Q,x-t',x'). 

In fact, gives us the particle mean number N{(3, /i) — 

As discussed in , the chemical potential fj, can also be included in the boundary condition for the Euclidean time 
direction, through the invertible transformation 

D{P, R) - ifj.j° = e'^'^^DiP, R) e-^'^*, 

(8) 

Khc{t, x; t', x') = e^^P^Ht, x; t', x') e'f''^*', 
where k{t,x]t' ,x') is the Green's function of the modified problem, 

D{f3, R) k{t, x; t', x') = 6{x - x')5{t - t') 

Be'"3*fc(t, x; t', x') = 0, for \x\ = 1, (9) 
fc(l, x; t', x') + e-i'l^kiO, x; t' , x') = 0. 

This function has the development 

oc 

k{t,x;t',x')=R J2 knix,x') e-'^-^^'-*'\ (10) 



with the frequencies given by 



(2n+ l)7r 



(11) 



, for n e Z, 



and kn(x]x') satisfying 



rj„i?ip^ (g) I2 - «) CT • V fc„(x;x') = 2;')^" ® I2, 



(12) 



Bkn{x;x') = 0. 

It is convenient to define dimensionless parameters r2„ = R^n, P- — RfJ- and z — R/p. 



(13) 



The solution of (^) can be expressed as the sum of a solution of the inhomogeneous equation, fco(i, x; t' , x'), with a 
vanishing behavior at infinity, and a solution of the homogeneous equation, k{t, x;t' , x'), chosen in such a way that 
the complete Green's function k{t, x; t' , x') obeys the boundary condition in (|^ ) . Th is sp littin g of the Green's function 
induces a similar separation for the mean fermionic number. In subsections II A| and [IB, we will considered both 
pieces separately. 



A. Fermionic mean number for the boundary less case 



In this case the mean fermionic number is given by 



No 



dG 



1, 



(14) 



where feg corresponds to the Green's function of the modified problem (^ , satisfying the temporal boundary condition 
and vanishing at spatial infinity. 



ko{t,x;t\x') ^ R J2 C (a;-^;') e 



-ifl^lSit-t') 



n— — 00 



k^^\x,x') 



1 e 



-s„n,Ax-x' 



An 



(15) 



In the previous expression, Sn — sign(n + 1/2). 

In order to get A'o we have to compute the trace in eq. (|lj). Due to the singular behavior of fco at the diagonal, 
we can take t' = t, but we should keep x' ^ x wp to the end of the calculation. In this way, 



No^j^j dt J r^dr / dQ { ti 



1 \ „ — SnQn\S:—x'\ 



'? o 



Z7 



— x'^ 



(16) 



The second term in the integrand of (|16|) vanishes because of the properties of Dirac matrices. For the first one, with 
e = a; — x' , we get 
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No = lim -^^ V a.e-^"^" 

n—~oc 



8 R d 



n=0 



So, 



47r / _ 2 M' 



(17) 



(18) 



in agreement with the value for a free fermionic gas, as can be seen, for instance, in Ref. p8[ . 

B. Influence of the boundary conditions on the fermionic mean number 

For the solution of the homogeneous equation required to adjust the boundary conditions for the full Green's 
function, we propose the development 



(19) 



fe=0 



and write A;q"''(x, x') as 



(x, x') =J2p''®A'' (x, x') . 

k=0 



(20) 



Concerning the A'^ coefhcients, A" and vanish, since they are the solutions of an everywhere regular homogeneous 
problem. For the other coefficients we have 



A^x,x')^-^A^x,x'), 

tl ill 



(21) 



+ {innY] A^{x,x') = innS'-^\x,x')u. 

Following |l|], we get A^ as an expansion in terms of spinorial spherical harmonics, 

A^{x,x')^ J2 4M^™(^'^')'/'J^^.™(^^)'/'^^^^™(^^'), 



(22) 



(23) 



where 



^Ui' mii^^'^') = 5i,i'iSn^'i3i{iSn^nr<)h']}\iSn^nr>) ■ 



In eq. (24), ji{x) and h\^'{x) are the spherical Bessel and Hankcl fmictions, respectively, and 



(r> = Max(r,r')). 

For the coefficients a'^ in eq. (|l9|), we get 



(24) 
Min(r, r') 



while the boundary conditions at r = i? imply 



= 0, 



\2" 



a" = 0, 



(25) 



a -fa} - 10? = iA^. 



Using the following Ansatz, 



we finally get 



j,l,l' ,m 



(26) 



(27) 



where Cj and dj are given by 



jj+l/2{iSn^n)hf_l-i^/2('^Sn^n) ^ jj-l/2iiSn^n)h'^''J_i/2{iSn^n) 
•^'1 + 1/2 (^'^" ") j 4 — 1/2^^ ^n^n) 



(1) 



(28) 



The correction to the fermionic mean number due to the boundary conditions is 



1, 



Ar=-Tr 



(29) 



(30) 



Once again, in order to get a well defined result, we must proceed carefully when treating the kernel of this operator 
at the diagonal. We can take t' = t and O' = O safely, while keeping r' = r(l — e), with e > 0. We thus get 



N 



00 „i 
2z lim iSn{i^n)^''^^Cj I 



drx 



jj+l/2{iSn^nr)jj+i/2{iS„^lnr')+jj_i/2{iS„^lnr)jj_i/2{iSn^nr')^^_^^^^ 

After a straightforward calculation of the radial integral, we obtain the following expresion 

00 

iV = lim Nn,k{z,p.;e) 

e— >0+ 



(31) 



00 'OC' 



= -2i.lim ^ ^ - J- (2. + 1) " + 'y-:!r X 



n=-ook=0 



{du,n 



(32) 
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where 



= 5„a.(l-e)^log/.(a;) 



ly = k + 



(33) 



a;=S„0„(l-£) 



and d^,„ = ^t.nL^Q- 



The double series ( p2D is not absolutely convergent for e = 0, so we must keep e > in the general term. A way to 
isolate such singular behavior consists in the substraction of the first terms in the asymptotic (Debye) expansion of 
Nk^n, ^k'n (^^^' fo"" example, The difference 



|iVfe,„(z,/i;e)-A*f„(z,/2;e)| 



„M+1 



(34) 



where = + can be considered as a "radial" variable in the double series. 
So, we can write 

TV = TVi + iV2 = 



oo oo 



(35) 



n= — oo fe=0 



ri= — oo fe=0 



In the first term of (p5|), for M > 3, the limit e ^ can be taken for the general term, inside the (absolutely and 
uniformly convergent jdouble sum. The nonregular behavior of the series is isolated in the second term of (|35|). 

Our calculation strategy consists in the numerical evaluation of the first term of (|3^), and the analytical study of 
the second one. Notice that we may take M > 3 in order to improve the numerical convergence of the first term. 



Asymptotic (Debye) expansions 
Introducing the Debye expansion for the modified Bessel functions , we can write 

with 



{-2iz)[\-e) 



(36) 



1 



/ ^+(^,^/(-l/2)(e,t)) 



t 



[S+{,y,t)S-{:.,t) - 5+(i.,i)5-(j.,i)+ t {S+{,.,t)S-{,y,t) - 5-(j.,i)5+(i.,t))] 



l-t2 

where stands for a consistent expansion up to the order /0~*^, and 



Eiu,t,e) = exp (^^ (/i/2(e,<) - l) + J/ log (^(1 - e) 



t+l 



/2(e,t) 



In the expressions above. 



M 



(37) 



(38) 



V 

t = - . 



/(6,t) = l-e(2-6)(l-t2 
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Ufe(t) 



oo 



.Ukit) 



k=l 

oo 



Vk{t) 



k=l 

oo 



fe=i 



(39) 
(40) 
(41) 
(42) 



with Uk (t) and Vk (t) being the polynomials appearing in the Debye expansion of Bessel functions, as defined as in 
Ref. lo). 

It is easy to see that 



i?(z/,i,e) = exp (-6^) (l + 0(e2)) 



n— — OO k—0 



No 



OO oo 



.n= — oo fe=0 



+0+ 



The first term, A^i, will be numerically evaluated from the expression 

Nn,k{z,fi,e = 0) = 



(43) 



So, the corrections to the fermionic mean number due to the boundary conditions (35) can be written as 

oo oo 

Ni= E [^".fc(^'^;^ = 0)-f^„D(j.,i,e = 0)] , (44) 



(45) 



(46) 



and from the asymptotic expansion 13, at e = 0. For computational convenience we take M — 6; the corresponding 
expression for D is 



D{i^,t,e = 0) ~2iz 



2 + 



"4 2" ^ ~ ^ T" 



■—r \ \ \ 3t^" 

z^4l8 416 8 2 



16 



1 f-t^ 3<6 27 1'' 483 1'^ 25*1° 1211 285*^2 781 213 

ll^ 8 32~ 32 ~8 32 16~ 32 16 



1 /-t*^ 47 2909^9 473 1^" 531 

1^ [i^ 64 r ^ 



256 16 



64 



7 



4255^12 14795^13 6297^1-* 481 1491 32799 



32 



128 32 



16 



256 



(47) 



The second term, N2, will be studied analytically in the following. Taking into acount (p3), we have 



00 00 



J2 ^n„i^(^^,t,e = 0)e-^''/* (1 + 0(6)) 



_n— — 00 k—0 



-i €^0 + 

Now, employing (E^), Equation (Esh can be naturally expressed in terms of the double series 



00 00 



n= — oo k=0 

(where, for convenience, we omit the dependence on fl and z) obtaining 

-s(5,0;e) 3s(5,l;e) 3s(5,2;e) s(5,3;e) 



N2 = lim 



16 



(48) 



(49) 



s(6,0;e) _ 3s(6,l;£) _ 3s(6,2;e) _ s(6,3;e) 
16 8 4 2 



47s(7, l;e) 27s(7,2;e) lls(7,3;e) 5s(7,4;e) 



64 



32 



16 



3s(8,2;e) 9s(8,3;e) 27s(8,4;e) 3s(8,5;e) 



2909s(9,3;e) 483s(9,4;e) 7s(9,5;e) 473s(10,4;e) 



256 



32 



16 



25 5(10. 5; e) „ „ , 531s(ll,5;e) 1211 s(ll, 6; e) 
3s(10, 6; e) — 



64 



32 



71s(ll,7;e) 4255 s(12, 6; e) 285s(12,7;e) 14795 s(13, 7; e) 



16 



32 



16 



128 



781s(13,8;e) 6297 s(14, 8; e) 213s(14,9;e) 481s(15,9;e) 



32 



32 



16 



1491s(16,10;e) 32799 s(17, 11; e) 



16 



256 



(50) 



® 



In the following section we will prove that s{p, q: e) has a finite e —t 0^-limit, which allows for the evaluation of 
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Evaluation of the double series s{p,q\e) 
In this section we study the double sums defined in equation (^9|), which can be written as 

oo 

where 



k=0 



It is easily seen that the Sp{h', e) are odd- functions of /I, and satisfy the recursion relations 

d 



(51) 



(52) 



(53) 



Sp{i^, e) 



1 



m- 2 



■ 1 a5p_2(i/,e) 

eSp-i[i^, e) + - 



v dv 



for m 7^ 2 , 



(54) 



which allow the computation of any Sp{v^ e) from the knowledge of S2{i', e)- The last series will be determined in what 
follows. 

Making use of the Poisson formula |^ , 52 {v, e) can be transformed into the series 



(55) 



where the coefficients are given by 



(-1) 



23 / dx ^_^^^e'^^/^e-'V(^^^^^^ 



(56) 



It is convenient to treat the terms with t — > Q and £ < separately. The corresponding integrals can be 
exactly computed in each case, by adequately extending the integration path on the complex x-plane. For /2 > 0, one 
gets 



Co = H{fx - v) 



1 - 



dx- 



— - sin {^\/ x^ — 



(57) 



Jy X-^ -U-^ J J 



(58) 



2 r°° T - 

C,<o = (-1)' <! -H{fi - v) / dx- e-'^P-^^'^ sin(eVa;2 - v"^) + 



TT /r. — 



iJ(i.-/i) 

where stands for the Hcavisidc step function, 



^-t{f.-u)/z ^ £_ / rfa; e-^(M-^)A sin(e V^^^^) 

TT ./., X-^ - 



(59) 
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Finally, the sum over the index i can be perfomed to get 

1 



1 



1 + e^(^^'')/^ 1 + e(^+'')/^ 



du 



TT ./n U 



sin(ei/u) . 



(60) 



Now, the expressions for 5p^2('^, e) can be inmediately obtained employing ( |53| ) and ( p4| 



All 5p(z^, e) so defined are regular at e = 0, and exponentially vanishing with v (even at e = 0). Such behavior 
guaranties the uniform convergence of the fc-series in expression (|5l|), allowing to take the e — > limit of Sp{v,e) 
therein, 



s(p,g;e = 0) = ^z.?5p(i.,e = 0). 



fc=0 



This limit also simplifies the recursion relations leading to 

1 



52k (1^,0) 



[2(k- 1)]!! \ vdv 



1 d 



(61) 



(62) 



where 



52.+i(^,0)= ^ {--^\ 5i(:.,0) 



[2k- 1]!! V vdv 



Si{v,Q) = — / du 

Jo 



52(^^,0) = 



1 + e^f^^*^)/^ 1 + e(^+'')/^ 



(63) 

(64) 
(65) 



Equations (|62|) to ( pSf) lead to exact expressions for iSp(z^, 0), allowing for the numerical evaluation of s{p, q;e = 0). 
We can also obtain analytically an expansion of s{p, q;e = 0) for large values of z = RT. In fact, for p > 4 it is 
straightforward to get 



s(p,(7;e = 0) = (-2z)-1^^25Ra, 



P-2,9 : 



where 



C<j oc 



fc=0 n=0 



Making use of the integral representation 



Pi 2 

dTT'^-^e-^P 



we obtain 



(66) 



(67) 



\fc=0 / \n=0 



(68) 



In this expression, for large z and for each n, high values of t are strongly suppressed by the exponential factor 
g-T(2irz)(n+2 ) ^ So, wc Can usc the Euler-Mac Laurin expansion [pol to evaluate the factor 
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E 



1 1 ^(d+i\ , 



2 ^ (^+ 1) 



-J - 2Bi+i 



- V ■ V - . - ^ (^+ 1)! 
fc=o " ^ ' e=i ^ ' 

Replacing in (^8|), solving the r-integral and the remaining n-sum, we get 



(69) 



(-2z)_a 

p — 2 dp. 



5ft 



/ 1 y-'^-3 r(£^)r(^) 



for p > 4. 



.^(-1)7 1 

ei \2ttz 



p-2+2e r ( p-^+^f^^ 



r(H^) 



C p-2 



III. NUMERICAL EVALUATIONS AND APPLICATIONS 



(70) 



The results obtained in the previous sections allow for the complete evaluation of the mean fermionic number and 
the /X— dependent part of the Gibbs free energy for the fermionic field we are considering. 

In order to obtain N, at fixed z, we can numerically evaluate Ni from equations (|4^), ( ]4^ ) and (^^, and N2 from 
(|^) and (je^) to (|6|). 

Adding these contributions to Nf), equation dl^), we obtain the fermionic mean number. The result is shown in 
Fig. for z = 1/8,1/4 



5 



1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

— z=l/8 
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/ 




- - z=l/4 


y 






// 






// 




























1 .... 1 .... 1 . 





0.5 1 1.5 2 2.5 3 

|LlR 

FIG. 1. Fermionic mean number for z = 1/8, 1/4 
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Notice that N{z, ^) has the expected behavior with /x, since it shows steps at the adimensionahzed eigenvahies of 
the Dirac Hamihonian ||2l[| . 

Numerical evaluations become more difficult when z grows up. But, in the z > 1 region, the asymptotic approxi- 
mation for N2 derived from ( |70| ) can be used. Moreover, in this region Ni decreases with z. In fact, the first term 

of (|70|) suggests that Ni contributes with a term of order z^"^. So, up to this order, we obtain, 



N 



fiz 



7r2/ Stt'^ 1267r3z2 



(71) 



The first two terms in the right hand side of ([7l|), proportional to and R respectively, reproduce equation (27) 
in Ref. where only the first finite size corrections to the free energy are taken into account. In the approximation 
leading to ( |7l| ) we were able to get the next, 1/R, term. 



Finally, the Gibbs free energy can be obtained by integrating the mean number. 



R (G(z, /2) - F{z)) = - / N{z, p,/) djii , 



(72) 



where F(z) stands for the Hemholtz free energy of the system. In fact, F{z) has been studied in Ref. g, so the 
expression above completely determines G{z,ji). 

If the fermionic number is fixed, in media, to N, one can obtain (numerically) the chemical potencial as a function 
of z, fi{z, N). Replacing in (f2|) we get 



RG{z,fl{z,N)) 



p.{z.N) 



dflfN{z,flf) - RF(z) 



(73) 



The numerical evaluation of ([7^ ) for N = 1/2, in a wide range of values of z, is shown in Fig. 




FIG. 2. Gibbs free energy at fixed fermionic number N — 1/2 
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Application to the model of quark-gluon droplet 



We will now apply our results to the study a simple model of hadrons: a static spherical droplet of confined quarks 
and giuons. 

Considering the degeneracies of flavor, Nf — 2, and color, Nc = 3, we obtain, from ( [73|) 



G 



quarks 



NfNc 
R 



fl{z,N=1,) 



d^JilN{z, Hi) + z(-/3Ftcrmions(z)) 



(74) 



where — /J-Ffcrmionsl^) is given in section III of Ref. 

On the other hand, in the 1-loop approximation, the gluon field can be described as 8 copies {Ng — 8) of the Abelian 
gauge field studied in Section IV of Ref. So, 



Co 



' R ^ ^ /^-^gaugc) 



(75) 



Finally, a bag stabilizing term, BV, is introduced. The bag constant is fixed to insure the pressure equilibrium 
when temperature reaches the critical value (which we take as Tq = ISOMeV). From the leading volume terms, we 
obtain 



90 ° 



(76) 



When the Casimir energies of quarks and giuons (obtained in Ref. in the framework of an analytic regularization), 
which can be written as C/ R with an undetermined coefficient C, are not considered, the remaining part of the total 
Gibbs free energy. 



G{Z, fL{z, N^3))^ Gquarks + Ggluons + BV 

is as shown in Fig. |3| for different values of the temperature. 

3000 
2000 
^ 1000 
S 
-1000 
-2000 



(77) 




0.5 



2.5 3 3.5 



R[fm] 



FIG. 3. Gibbs free energy for the quark-gluon droplet. From top to bottom: T = 140, 150, 160, 170MeV 
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The situation we find is consistent with the first order transition of a hadronic bag to a deconfined quark-gluon 
phase described in For each positive value of C, and for < T < Tq, an absolute minimum of G{z,jl{z^N = 3)) 
appears at a finite value of R, corresponding to a stable bag configuration. When the temperature reaches Tqi this 
becomes only a local minimum representing a metastable state. When T grows up, the local minimum and the local 
maximum approach each other. Finally, there is a critical temperature Tc > Tq such that for T > Tc there is no local 
minima at finite i?, and only the deconfined phase can exist. 

The computational techniques employed in this paper can, in principle, also be applied to more refined models of 
hadrons, such as the hybrid chiral bag model [ p2[p^ . In Ref. the behavior of the total energy at T = and 
fi ~ has been analyzed, showing a remarkable independence on the bag radius, in agreement with the so called 
Cheshire Cat Principle This justifies the proposal in [|l5|,^, of looking for a deconfining transition only at the 
temperature (and ^) dependent part of the free energy. 

The Helmholtz free energy (/i = 0) of a chiral bag can be obtained from the results in |^,^. The method developed 
in the present paper could be applied to study the ^ dependent part of the Gibbs free energy for this system. Notice 
that, for this model, boundary conditions (which depend on the chiral angle) make the evaluation of the Green's 
function more difficult, since there symmetry corresponds to the diagonal subgroup of 5'[/(2)i.otationai ® 'S'?/(2)isospin 
(see Ref. for the T = case). We will report on this subject elsewhere. 
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